Sublocalization, Superlocalization, and 
Violation of Standard Single Parameter Scaling in the Anderson Model 



Jan W. Kantelhardt^^^ and Armin Bunde^ 

^ Institut fUr Theoretische Physik III, Justus-Liebig-Universitdt Giessen, D-35392 Giessen, Germany 
^ Center for Polymer Studies and Department of Physics, Boston University, Boston, MA 02215, USA 
(submitted: January 19, 2002, revised version: April 24, 2002) 

We discuss the localization behavior of localized electronic wave functions in the one- and two-dimensional 
tight-binding Anderson model with diagonal disorder. We find that the distributions of the local wave function 
amplitudes at fixed distances from the localization center are well approximated by log-normal fits which be- 
come exact at large distances. These fits are consistent with the standard single parameter scaling theory for the 
Anderson model in Id, but they suggest that a second parameter is required to describe the scaling behavior of 
the amplitude fluctuations in 2d. From the log-normal distributions we calculate analytically the decay of the 
mean wave functions. For short distances from the localization center we find stretched exponential localization 
("sublocalization") in both, d = 1 and d = 2. In d = 1, for large distances, the mean wave functions depend 
on the number of configurations TV used in the averaging procedure and decay faster that exponentially ("super- 
localization") converging to simple exponential behavior only in the asymptotic limit. In d = 2, in contrast, the 
localization length increases logarithmically with the distance from the localization center and sublocalization 
occurs also in the second regime. The A'^-dependence of the mean wave functions is weak. The analytical result 
agrees remarkably well with the numerical calculations. 

PACS numbers: 71.23.An, 73.20.Fz 



I. INTRODUCTION 

The Anderson model is the standaijd model for the elec- 
tronic, properties of disordered solids,El for reviews, see e. g. 
Refs.l3a The wave functions are known to be localized in di- 
mensions d < 2 (except at the band center for non-diagonal 
disorder) as well as for sufficiently strong disorder in d > 2. 
The shape of these localized wave functions is believed to be 
characterized by an asymptotic alb' exponential decay, which 
can be described most effectivelyB by the Lyapunov exponent 
7 or its inverse, the localization length A. The simple exponen- 
tial shape has been the basis of many calculations and models, 
e. g. Mott's hopping conductivity and tunneling on quantum 

, dots. Even for other models of disordered matter, e. g. for per- 

' eolation models,! simple exponential localization of the wave 
functions has been assumed based on the analogy with the 
Anderson model. 

In addition to the exponential decay at large distances from 
the localization center, the wave functions in the Anderson 

' model exhibit large fluctuations of their amplitudes.B These 
fluctuations are rather important for transport properties, e. g. 
for the distribution of mesoscopic conductances, since the 

, core region of the wave functions is dominated by the fluc- 
tuations, and the exponential decay can only be observed at 
distances much larger than the localization length A for situ 
gle wave functions. In the single-parameter scaling theorydiJ 
(SPST) it has been assumed that the whole distribution of the 
conductances or Lyapunov exponents for an ensemble of con- 
figurations can be described by one parameter only, e. g. the 
width of the distribution of values of Lyapunov exponents de- 
pends on the average Lyapunov exponent 7 in a universal way. 
Hence, the localization length A = I/7 could be used as sin- 
gle relevant scaling parameter fully characterizing the local- 
ization behavior. 



Very recently, there has been a growing interest in the ques- 
tion, if the behavior of the fluctuations of the wave funciions 
are nnivprsal apH can be described within the SPSTJjiIj see 
also Refs.ElOiij The interest is partly motivated by the ex- 
perimental observations of a transition in the behavior of the 
conductance of low density two dimensional semiconductor 
devices from an insulating like temperature dependence at low 
densities to a metallic one at higher densities.E^lThis apparent 
localization-delocalization transition in d = 2 (2d) is at odd 
with the SPST for non-interacting electrons. 

In this paper, we investigate the distribution of the ampli- 
tudes of Anderson locaUzed wave function in d = 1 (Id) and 
2d at given distances r from the localization center, which al- 
lows a test of the SPST hypothesis. We show that the SPST 
hypothesis holds perfectly in Id, while deviations from the 
standard SPST occur for localized wave functions for the An- 
derson model in 2d. Deviations also occur for the percolation 
model in 2d and 3d. In both cases, the width of the ampli- 
tude distribution is characterized by a second parameter that 
depends on the type of disorder In addition, we find a clear 
logarithmic system-size dependence of the localization length 
A for the Anderson model in 2d, but no significant dependence 
on the boundary conditions. 

From the amplitude distributions we can calculate analyti- 
cally the behavior of the mean wave functions at sufficiently 
large distances from the localization center. We find that sim- 
ple exponential decay as usually anticipated occurs only in Id 
in the asymptotic limit where the wave function has decayed 
to extremely small values. At finite distances from the local- 
ization center, we observe two different localization regimes. 
In the first regime, at short distances from the localization cen- 
ter, we find stretched exponential localization with a localiza- 
tion exponent below one ("sublocalization") in Id and 2d. In 
the second regime, at larger distances, the averages in Id de- 
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pend on the number of configurations and we observe "super- 
localization", where the localization exponent is greater than 
one. The crossover from the sublocalization to the superlocal- 
ization regime increases logarithmically with the number of 
configurations. In the asymptotic limit, the localization ex- 
ponent converges to one, yielding the theoretically predicted 
simple exponential behavior In c? — 2, however, we observe 
the surprising result that the sublocalization of the mean wave 
functions remains also in the second regime. 

The paper is organized as follows: In Section II we briefly 
introduce the Anderson model and summarize some of the 
known results. In Section III we discuss distributions of wave 
function amplitudes, beginning with a review of the predic- 
tions of the standard SPST and the random matrix theory. 
Then we introduce an amplitude statistics which is most use- 
ful for the characterization of localized wave functions, and 
present our numerical results for this amplitude statistics. In 
Section IV we calculate analytically the decay of the mean 
wave function amplitudes and distinguish different localiza- 
tion regimes. After describing an appropriate averaging pro- 
cedure, we present our numerical results for the decay of the 
wave functions in Section V. We summarize our results in Sec- 
tion VI. 



II. THE ANDERSON MODEL 

We consider the Schrodinger equation in tight-binding ap- 
proximation for the wave function of a quantum particle on 
regular lattices. The coefficients for each lattice site n sat- 
isfy the tight-binding equation 

Elpn ^ eni^n +^^Vn,n+5ll^n+S, (1) 
5 

where the sum runs over all nearest neighbor sites n + S of 
site n. IV'nP is the probability density on site n for the wave 
function with energy E. The hopping terms Vn^n+s are con- 
stant for nearest neighbor cluster sites and zero otherwise; for 
simplicity we take V ^ 1 as energy unit. The values for 
the diagonal terms e„ are randomly chosen from the inter- 
val [— w/2, w/2], where w is the disorder parameter (diagonal 
disorder). 

Here, we consider the Anderson model with periodic as 
well as hard wall boundary conditions for a linear chain and on 
the square lattice.^ Iti^ese cases the wave functions are known 
to be localizedou LrllJ even for small disorder strength w. For 
large distances r from the localization center, the decay of the 
eigenfunction amplitudes is believed to be exponential. Thus, 
a localization length A and the Lyapunov exponent 7 = A^^ 
can be defined to describe the asymptotic decay of the mean 
logarithm of the amplitudes, 

(ln|V;(r)/V(0)|> =-r/A = -r7 (2) 

for r ^ X. In Id the localization length in the band center 
i? ~ can be calculated analytically from a recutsion equa- 
tion derived from Eq. ([l|) using Greens functions:ta A(w) ~ 



105.045/w^. This result has also been verified numerically.ta 
While the prefactor slightly depends on the energy E, the 
l/w^-dependence holds for all eigenstates except at-the band 
edges, where a different analytical theory appliesBO In 2d 
there is no analytical formula for the localization length A. 
The numerical results published so far were obtained using 
finite size scaling combined with (a) transfer matrix calcula- 
tions of the cojoductancetJ and (b) exact diagonalization and 
level statistics.Eil 

We are interested in the distribution of the amplitudes |i/)„| 
for the eigenfunctions of Eq. Because of the disorder, 
the eigenvalue equation has to, be solved numerically. In Id 
we use aB-Jteration method,E3 while we employ a Lanczos 
algorithmcd in 2d. 

III. THE AMPLITUDE DISTRIBUTIONS 
A. Definition 

In the original single-parameter scaling theoryB (SPST) 
the fluctuations of the wave functions had not been consid- 
ered. An exiension of the theory includes the scaling of the 
fluctuations. □ Still it has been assumed that there is only one 
relevant scaling parameter characterizing the localization be- 
havior and the fluctuation behavior. This means that, e. g., the 
position of the center as well as the width of the distributions 
of conductances, Lyapunov exponents, or localization lengths 
can be characterized by only one scaling variable. Practically 
the distribution of the parameter 

7(i) = ^ln(l + l/5) (3) 

with the conductance g and the system size L has been inves- 
tigated theoreticaJJyO within the so-called random-phase hy- 
pothesis, see Ref.B for aJhorough discussion. While the aver- 
age (7(i)) was shownU to approach the Lyapunov exponent 

7. 

(7(L))-7=A-i for L » 7"', (4) 
its fluctuations were expected to become 

(7'(L)) - (7(L))2 = 7/L for L » 7-1 (5) 

asymptotically. Thus, 7 fully characterizes the asymptotic 
scaling behavior of the average j{L) and its fluctuations. 
Although many numerical calculations of (jiL)) [testing 
Eq. (^] have been reported, there is only little work devoted 
to the scaling behavior of the fluctuations and their relation to 
7. Recently Deych et alJS considered the Lloyd model (which 
is identical to the Anderson model with diagonal disorder ex- 
cept for a Cauchy distribution of the site energies) and found 
that states at the tail of the density of states violate the SPST 
and that there exists a new length scale that is responsible for 
this violation. 

In this paper, we determine numerically the distribu- 
tion of the logarithm of normalized amplitudes, A{r) = 
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— In |'0(r)/-0(O)|, for large distances r from the localization 
center (at r = 0). The distance r is directly related to the sys- 
tem size L used in Eqs. (||)-(^, since the value and the fluctu- 
ations of the conductance g for a system of size L reflect the 
value and the fluctuations of its wave functions at the distance 
r = L from the localization center. Following Eqs. and 
(0), it is obvious that 



(A(r))^(-ln 



= _ = r(7(r)) 



(6) 



for r 3> A. Hence, A{r) corresponds to rj{r), and the fluctu- 
ations of A{r) are expected to become 

{A\r)) - {A{r)f = ((f) - (7)^)^^ = 7^ = r/X (7) 

for large r in the standard SPST with the random-phase hy- 
pothesis following Eq. In order to estimate possible de- 
viations from this standard single parameter scaling behavior, 
we will calculate the fluctuation parameter 



a{r)^2[{A^r))-{A{r)r] A(r)A 



(8) 



which is expected to become a = 2 asymptotically in the stan- 
dard SPSTeI according to Ea. for the Anderson model and 
(7 = 4 for the Lloyd model .B 

In an alternative approach it hasJpeen shown some years 
ago, that the random matrix theoryo can be applied_to elec- 
tronic eigenfunctions in disordered systems in 3dE3 and in 
2d!E3 In 2d, for small disorder strength w, the localization 
length becomes larger than the system size, A 3> -L, and 
the wave functions become quasi extended. In this regime, 
as well as in the metallic regime in 3d, the random matrix 
theory correctly describes the level spacing distribution as 
well as the distribution of the eigenfunction amplitudes | ipn \ , 
which are usually assumed to be normalized according to 
X)n IV'nP = l^p^ y ipt ^ity distribution (histogram) is usu- 
ally defined as 



fit) 



(9) 



to be consistent with the numerical results for real fluctuating 
Anderson wave functions.Eil Thus, a simple exponential de- 
cay without any fluctuations cannot be distinguished from a 
real fluctuating wave function, if ^ is used to calculate the 
intensity statistics. The fluctuations simply do not affect f{t). 
This indicates that the definition (^ is not appropriate to fully 
characterize localized eigenfunctions. Further, since no refer- 
ence to the localization center is taken into account in (^), the 
calculation of the decay of mean wave functions is not possi- 
ble from /(t). Another problem is that f{t) depends on the 
system size N for localized states. So, we introduce a differ- 
ent definition for the amplitude distribution here, which avoids 
these problems. 

Since we are interested in the decay of the localized wave 
functions, we investigate the distribution of the amplitudes 
at a certain distance r from the localization center By defi- 
nition, the localization center for each eigenfunction {ipn ) 
is that site tiq, where the maximum of the amplitudes oc- 
curs, i. e. IV'no | = max„ Itpn \- An alternative definition 
for the localization center would be to consider the lattice 
point most close to the "center of gravity" of the wave func- 
tion as the center of localization. That is the site ni where 
J^nii^n - Xni)'^ + (Vn - yn^Y\\i>nV rcachcs its minimum 
value (xn and j/„ designate the coordinates of site n here). Of- 
ten, these two definitions of the localization center coincide. 
But even if they do not coincide, there is practically no dif- 
ference in the results for the amplitude distributions, as we 
will show below. In order to achieve a system-size indepen- 
dent normalization of the amplitudes which allows us to 
compare numerical data for different system sizes and simpli- 
fies the analytical description significantly, we choose to set 
If/jngl = 1. Note that this choice is different from the usual 
normalization IV'nP — 1 underlying the definition (^, 
since the sum is always larger than one for our normalization. 

To describe the actual decay of the localized wave func- 
tions, we do not consider the distribution of the amplitudes for 
the whole eigenfunction [as in the definition (^], but instead 
look at the distributions of the amplitudes at given distances 
r from the localization center Hence, we replace /(t) by the 
histogram distribution function (v4, r) defined by 



where the sum runs over all lattice sites, N is the total num- 
ber of lattice sites, and the average (•••)£; is taken over sev- 
eral eigenfunctions with energy E' w E. In the limit 
of small disorder strength m. f{t) is given by the universal 
Porter-Thomas distribution,EZI 



/pT(i) = 



:Cxp(-i/2). 



(10) 



This result has been confirmed numerically recently for quasi 
extended states in the-^jXaderson model in 2dt3 and for the 
metallic regime in 3d.ll3il3 

The intensity distribution f{t), however, cannot appropri- 
ately describe the behavior of localized wave functions in the 
limit of strong localization A <C i. The analytical form of 
/(t) obtained by inserting a simple exponential decay without 
any amplitude fluctuations into the definition (m was shown 



ff(Ar) = {^f]<5(r-r„)<5(A + ln|^(^')|)) , (11) 



where r„ is the distance of the site n from the localization 
center, and A^,. is the number of sites with r„ = r, i. e. 
Nr = J2n=i ^(^ ~ ^")- definition (pi]), we use the 

variable A = - In | V'i^'^ | instead of t = iV'i^'^ ^ from ^ 
for two reasons: (i) t is system-size dependent for localized 
states, and (ii) A can be related much more easily to the pre- 
dictions of the standard SPST, that we want to test. In fact, we 
have already derived the relation in Eqs. (^) and (^. Note 
that H{— \n\i/j\,r) is the probability density of finding the 
amplitude at distance r from the localization center In 
the following, we will use this definition of the eigenfunction 
amplitude distributions to present our numerical results and 
discuss the localization behavior. 
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FIG. 1.: Histogram distributions H{A,r) of the absolute ampli- 
tude values \tp\ at distance r from the localization center versus 
A = — In I?/; I for electronic wave functions in the Anderson model 
in Id with disorder strength w = 1.0 at several fixed distances r 
from the center of localization ranging from r = 1600 (leftmost 
curve) with increment Ar = 1600 to r = 28800 (rightmost curve). 
To determine the distributions, more than 10^ eigenfunctions have 
been calculated on chains of lengths L = 30000. 



B. Results in Id 



Figures |i] and ^ show our numerical results of the ampli- 
tude distribution functions H{~\n\ip\,r) for the Anderson 
model in Id with periodic boundary conditions. Using an iter- 
ation methodE3 we have calculated eigenfunctions with en- 
ergy E K (coiTesponding to the band center) for linear 
chains with five different disorder strengths w and determined 
the distribution functions H{A, r) for several values of r. Fig- 
ure |l] shows 18 of these histogram distributions calculated for 
more than 10^ chains of length L — 30000 with w = 1.0. 

We find that the histograms In j-i/'j, r) of the ampU- 
tudes j-^l at fixed distances r from the localization center sur- 
prisingly well follow a log-normal ansatz, 



iJi„(-ln|V|,r) = 



A 



exp 



{ln\ij\+r/Xy 



crr/X 



(12) 



In general, both parameters, A and cr, might depend on r. 
Since H\n{A,r) is a simple Gaussian distribution with aver- 
age value {A{r)) = J AH\n{A,r)dA = r/A and variance 
{A^ (r)) — {A{r))'^ = (Tr/2\, the two parameters A and a cor- 
respond to the localization length and the fluctuation parame- 
ter, respectively, according to Eqs. (||) and (^. Note that the 
amplitude | V'l = exp(— r/A) occurs with the largest probabil- 
ity at distance r from the center of localization, which further 
prompts us to identify the parameter A with the localization 
length. In the following, A(r) always refers to the parameter 
determined from the distribution H{A, r) for a given r, while 
A represents the asymptotic value of A(r). 
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FIG. 2.; Scaled distributions ffs(^s,r) = H{A,T)yJ nar / \ of 
the amplitude values IV^I at distance r from the localization center 
versus As = {A ~ r/\)/ ar jX for electronic wave functions in 
the Anderson model in Id with disorder strength w = 4.0, 2.0, 1.0, 
0.5, and 0.35 at several distances r from the center of localization 
ranging from r ~ 20A to 250A. Due to the scaling the data collapse 
onto single curves. The curves for different w have been shifted by 
multiples of 0.2. The continuous line below the data corresponds to 
the ansatz Eq. ([l2|). To determine the distributions, (a) 4,200,000, 
(b) 3,400,000, (c) 1,100,000, (d) 200,000 and (e) 170,000 eigen- 
functions have been calculated on chains of lengths L > 250A(it;). 

The log-normal distribution of the amplitudes can be mo- 
tivated by the following argument. If Id wave functions are 
calculated using the transfer matrix method, the value of 
for large n is related to the values of for the two initial 
sites, n = and n = 1, by the product of n transfer matrixes 



n 



E 



1 



(13) 



This product of matrices will involve many products of the 
random site energies e„. Products of many independently dis- 
tributed random numbers asymptotically follow a log-normal 
distribution according to the central limit theorem (under 
some conditions, which we do not want to discuss in detail 
here). Thus, it is reasonable that i\}n will also follow a log- 
normal distribution asymptotically for large n. 
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The agreement between the numerical distributions 
H{A,r) (shown in Fig. and the log-normal distribution iJjn 
[Eq. (|l2|)] becomes perfect for large r values. This indicates 
that the distribution H{A, r) follows Eq. dl^ ) asymptotically. 
In addition, we find that both parameters, A and a, are inde- 
pendent of r. Fits of Eq. (|l2|) to the numerical distributions 
H{A, r) for w = 1.0 shown in Fig. |]give A = 105 ± 1 and 
(7 = 2.1 ± 0.15. Tha value of A is in perfect agreement with 
the analytical resultllS \{w) ~ 105.045/w^. The value of a is 
consistent with a — 2 expected in the standard SPST 

Since both parameters do not depend on r, the distributions 
shown in Fig. ^can be rescaled in order to obtain a data col- 
lapse for the different distances r from the localization cen- 
ter. If we multiply H{A, r) by y^rtrr/A and replace A by 
As = (A — r / X) / y/arJX, all distribution will collapse upon 
each other This is shown in Fig. ||, where data for the other 
disorder strengths w we considered are also included. Since 
the dependence of the localization length A on w is known, 
and the fluctuation parameter a turns out to be independent of 
w, we are able to rescale the distributions for all r and w in 
order to obtain a full data collapse also for the other disorder 
strengths. The numerical result ct = 2.1 ± 0.15, which is con- 
sistent with a ~ 2, fully confirms the standard SPST for the 
Anderson model in Id. 



C. Results in 2d 

For the Anderson model in 2d the scaling behavior of the 
histograms of amplitudes at fixed (Euclidian) distance r 
from the localization center turns outJij be not as simple as 
in Id. Using the Lanczos algorithrro we have calculated 
eigenfunctions of Eq. (|l]) for energies E ^ and two dis- 
order strengths w = 8.5 and w — 10.0. For these values 
of w the amplitudes of the localized states decay roughly to 
10~^^ and 10^^", respectively, on our 300 x 300 square lat- 
tice. This allows to observe the localization effects over sev- 
eral orders of magnitude. Quadruplicate precision is used in 
all calculations in order to get reliable values also for small 
amplitudes \tlj\ « lO^^*'. We compare the results for peri- 
odic boundary conditions (PBC) and hard wall boundary con- 
ditions (HWBC). 

We find that the histogram distributions H{A, r) can still be 
well approximated by the log-normal distribution H\n{A,r) 
[Eq. (p^)]. But now the parameter A(r) that corresponds to 
the localization length clearly depends on r. This dependence 
is shown in the main part of Fig. |[ where A is plotted versus 
r on a log-linear scale for systems with the disorder strength 
w = 10.0 and PBC as well as for w = 8.5 and HWBC. We 
also checked w = 8.5 with PBC and w = 10.0 with HWBC 
and found no significant alterations in the behavior Further, 
the results do not change if we disregard all those eigenfunc- 
tions that did not reach a preset accuracy limit in the Lanczos 
calculations or select only eigenfunctions within a reduced en- 
ergy interval. The results also do not depend on the way the 
center of localization is defined (as maximum or "center of 
gravity", see Section III. A). 
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FIG. 3.: Localization length A(r) of Anderson wave functions in 
2d in the locahzed regime for (a) w — 10.0 and (b) w — 8.5. 
The values of A(r) have been determined by fits of Eq. (|l2|) to the 
amplitude distributions H{A, r) for each distance r from the local- 
ization center. Data for periodic boundary conditions (PBC, a) as 
well as hard-wall boundary conditions (HWBC, b) have been an- 
alyzed. The results for both definitions of the localization center 
[as the maximum (MAX) and as the "center of gravity" (CG)] are 
practically identical (except for r < 10) and well described by a 
logarithmic dependence of A on r: (a) A(r)/0.89 — ln(r/1.2) and 
(b) A(r)/1.56 = ln(r/2.6). The straight lines correspond to these 
fits, but they are shifted by a constant offset for visibility. In the 
inset, the estimated fluctuation parameters a{r) of the distributions 
are shown. They are approximately independent of r and w. We 
find a — 1.3 ± 0.2. For both disorder strengths, about 1000 eigen- 
functions were calculated with the Lanczos algorithm with quadru- 
plicate precision on a 300 x 300 lattice. 

The values of A(r) and a{r) shown in Fig. ||have been ob- 
tained by fitting Eq. ( p^ to H{A, r) for each r. Finite-size 
effects are clearly observable for PBC, while they seem to be 
rather negligible for HWBC. The rather strong finite-size ef- 
fects in the case of PBC are probably due to the many addi- 
tional paths leading from the localization center to a remote 
site via the edges of the lattice. In all cases the A(r) depen- 
dence is clearly a logarithmic one for r > 10: 



\{r)/cx = ln(r/cr), 



(14) 



where the two parameters cx and describe the scales of the 
localization length A and the distance r, respectively. Only for 
small r values (i. e. r <C 20), a power-law fit is also possible. 
The logarithmic dependence is observed over more than one 
order of magnitude in r in all our numerical simulations. 
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value clearly deviates from a = 2 expected for the standard 
SPST. 



TABLE 1 : The values of c\ and Cr from fits of Eq. ([l4|) to the local- 
ization length A(r) as well as the distribution width cr for the An- 
derson model in 2d. The results for two disorder strengths u; = 8.5 
and w = 10.0, for both types of boundary conditions (BC, HWBC 
= hard-wall, PBC = periodic boundary conditions), and for different 
subsets of the eigenfunctions are reported. The value of Smax de- 
termines the width of the energy interval [0, -Bmax] for the selected 
eigenfunctions. The number of eigenfunctions involved in each fit 
is also reported. The fitting range is 30 > r > 120 in all cases, and 
the maximum has been used to define the localization center. 

Table 1 reports the fitted values of cx and Cr for both disor- 
der strengths, w = 8.5 and w = 10.0, for both types of bound- 
ary conditions, and for different subsets of the eigenfunctions. 
The value of i?niax determines the width of the energy inter- 
val [0, E'max] for the selected eigenfunctions. The table shows 
that c\ and clearly depend on w, but not on the width of 
the energy interval. The values for different boundary condi- 
tions are also very close to each other, and we cannot find any 
systematic dependence upon the boundary conditions. Both, 
Cx and Cr, are rather close to the microscopic distance be- 
tween the sites (set to one here), but they are definitely not 
simply identical to the microscopic distance, since both de- 
pend strongly on the disorder parameter w. The significant 
logarithmic dependence of the localization length A(r) on the 
distance r from the localization center may explain the some- 
what different numerical resuJtsJor the localization length that 
have been obtained before.OO When comparing with those 
results one has to keep in mind, though, that they were deter- 
mined using finite size scaling of the conductance and level 
statistics, respectively, while our determination does not in- 
volve finite size scaling. In the finite size scaling techniques, 
one value of the localization length or a similar level statis- 
tics parameter is calculated for each configuration. Then the 
average finite size parameters are extrapolated to infinite sys- 
tem size by rescaling the results for different disorder w in a 
scaling plot, assuming the existence and scaling behavior of 
Aoo- Here, in contrast, we do not assume any scaling proper- 
ties of the localization length and determine the dependence 
of A(r) upon the distance r from the origin for each system. 

In contrast to A(r), the parameter a characterizing the width 
of the amplitude distributions and shown in the inset of Fig. || 
is still independent of r, except for the pronounced finite-size 
effects in the case of PBC. We find cr = 1.3 ± 0.2 for both 
disorder strengths w, both types of boundary conditions, and 
all selections of eigenfunctions as reported in Table 1. This 
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FIG. 4.: Scaled distributions Hs{As,r) = H{A,r)^J%ar/X{r) 
of the amplitude values at (Euclidian) distance r from the local- 
ization center versus yls = [A — r/\(r)\/ ^ ar/X{r) for electronic 
wave functions in the Anderson model in 2d with disorder strength 
(a) w = 10.0 (PBC) and (b) w = 8.5 (HWBC) at several fixed dis- 
tances r from the center of localization, (a) r = 8, . . . , 106, and (b) 
r = 20, . . . , 182. The r-dependent localization lengths \{r) for the 
rescaling procedure have been taken from logarithmic fits in Fig. 3, 
and the width parameter a is fixed to 1.35. For each part of the 
figure, ~ 1000 eigenfunctions were calculated using the Lanczos 
algorithm with quadruplicate precision on a 300 x 300 lattice. 

Using the logarithmic A(r) dependence, Eq. (0), with the 
parameters from Table 1 the amplitude distribution functions 
H{A, r) for the Anderson model in 2d can be rescaled to ob- 
tain a data collapse similar to Fig. ^ The result is shown 
in Fig. ^ again for w = 10.0 with PBC and w = 8.5 with 
HWBC. The numerical results confirm the log-normal shape 
of the amplitude distribution H{— \n\ip\,r). 

In summary, we find two major deviations from the stan- 
dard SPST for the Anderson model in 2d: (i) The logarithmic 
A(r) dependence indicates that an asymptotic value A = 
(assumed in the SPST) might be ill defined, since the loga- 
rithmic function does not converge, (ii) The fluctuations still 
follow the SPST hypothesis Eq. (^, if A is replaced by the 
r-dependent localization length A(r) in Eq. (|4|), and an ad- 
ditional prefactor cr/2 is inserted. But since this prefactor is 
smaller than one, the fluctuations in the Anderson model in 
2d are smaller than those predicted by the standard SPST hy- 
pothesis. 

In order to find out, if these two types of deviations from the 
standard SPST are a generic feature of any system in 2d, we 
have determined the amplitude distributions and the values of 
A and a for electronic eigenfunctions on percolation clusters 
at the critical concentration, which we had studied in a previ- 
ous work.L3 PercolationO is a standard model for structurally 
disordered solids. For site percolation, the Schrodinger equa- 
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tion in tight-binding approximation is identical to Eq. with 
e„ = for occupied sites (concentration p), and e„ = oo for 
unoccupied sites (concentration 1 — p).E3 The electronic wave 
functions are known to be localized for percolation clusters 
at the critical concentration p = pc for any embedding di- 
mension. We found that the amplitude distributions at given 
topological distances from the localization center again obey 
the log-normal distribution ( |l2[ ) for clusters in 2d as well as 
in 3d. A data collapse is obtained with constant parameter 
A, and no logarithmic dependence of A on the distance is ob- 
served. Thus, the first deviation from the standard SPST that 
we observed in the 2d Anderson model does not show up for 
the percolation model. Nevertheless, the second deviation, the 
modified value of the fluctuation parameter cr, also occurs in 
percolation. But in contrast to the reduced value cr = 1.3 for 
the Anderson model in 2d (compared with the standard SPST 
value a — 2), we find increased values w 2.6 and 3.0 for 
percolation in 2d and 3d, respectively, which are still smaller 
than (7 = 4 for the Lloyd model. We conclude that a model de- 
pendent second parameter a is required to describe the scaling 
behavior of the amplitude fluctuations. 



IV. ANALYTICAL RESULTS FOR THE MEAN WAVE 
FUNCTION 



are small and A = — \n jV'l large. Thus, we insert H\a{A, r) 
for H{A, r) in Eq. ([l6|). However, the exact shape of H{A, r) 
for small distances, corresponding to large amplitudes, is not 
known exactly. Hence, our result will be exact only for large 
values of r. We further have to assume that the number of 
configurations N is finite, so that the extreme tails of the dis- 
tribution H{A, r), i. e. the very rare events possibly involving 
small A, are not important. Our treatment includes a descrip- 
tion of the typical average over one configuration (one eigen- 
function), which corresponds to = 1. 

For a finite number N of configurations, the total number 
of sites at distance r from the localization center is identical 
to NNr with Nr = ar'^~^ and a = 2 (27r) for the Ander- 
son model in Id (2d). Clearly, those values of \'4>\ with a too 
small probability, i. e. probability smaller than 1/(A^ Nr), 
are unlikely to occur in N typical configurations. Thus, we 
must limit the distribution H\n{A^r) by cutting off the area 
1/{2N Nr ) at each of the two tails. According to Eq. the 
area below a lower cutoff A^iin is given by 



Hi„(A,r) dA 



1 



1 + erf 



(17) 



Since we can very well approximate the distribution of 
eigenfunction amplitudes \ip\ at large distances r ^ A by the 
log-normal distribution ([l2|), we can directly calculate the de- 
cay of the average wave function for large r by sim- 
ple integration. We can also estimate the way this average, 
which we denote by ^jv(r), depends on the number of aver- 
aged configurations N . In particular, we can obtain Nir) 
in the asymptotic regime which is not accessible numerically. 
The calculation is valid for Anderson wave functions in Id and 
2d (as well as for the percolation model), since the same log- 
normal distribution Hia{A, r), applies to all cases and only the 
two parameters A and a have to be adapted. Note, that the r- 
dependent localization length A(r) from Fig. ^has to be used 
for the Anderson model in 2d, while A is constant otherwise. 

In order to describe the localization behavior properly, we 
are led to the ansatz 



where erf (x) = 27r JJ^ exp(— i^) dt is the error function. 
Equating this area to 1/(2A^ Nr) = {2Nar'^^^)^^ and using 
the inverse error function erfinv(a;) we obtain the lower cutoff 
value 



A^ 



max < 



— erfinv f 1 -t—t 

A V Nar<^-^ 



(18) 



(15) 



which replaces the lower integration bound in Eq. (p^. 
Hence, for finite N, Eq. (|l6|) becomes 

/>oo 

^N{r)^ eM-A)H,^{A,r)dA. (19) 

"'A„i„(r,Af) 

The integration can be performed straightforwardly and gives 



where d^ is an effective localization exponent. If we average 
over all configurations, the resulting quantity ^oo is related 
to the distribution H{A, r) by 



*oc(r) = mr)\) 



exp(-A)iJ(A,r) (16) 



with A = — In I?/) I . Note that only positive values of A are 
possible because of our normalization ?/;„q = 1 for the local- 
ization center (maximum) at site riQ. In the preceeding Section 
we have numerically shown that H{A, r) can be well approx- 
imated by the log-normal distribution H\n{A, r) for large dis- 
tances r from the localization center, where the amplitudes 



'^N[r) = i exp 



- 1- 



4/ A 



1 — erf 



erfinv | 1 -r—r 

4A V Nar<^-^ 



(20) 



The first term of this equation seems to indicate a regular ex- 
ponential decay, if the second term could be disregarded and A 
was constant. But, indeed, the second term, involving the er- 
ror function as well as a complicated dependence on the num- 
ber of configurations N, cannot be disregarded. Thus, the 
average wave function '^N{r) depends explicitly on N and 
exhibits significant deviations from the exponential from. 
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FIG. 5.: The theoretical effective localization exponents d<s, de- 
termined by fits of the ansatz ( p^ to the analytical result ( po[ ) are 
shown versus r/X for several numbers of averaged configurations 
A'^ (see legend). For the Id results (upper curves) the parameters 
A = 105 and cr = 2 have been used, while A(r)/0.89 = ln(r/1.2) 
(from Fig. ^ and a = 1.3 for the 2d case. In the Id case we find 
strong dependence on A'^ and > 1, indicating superlocalization, 
while dxi < 1 indicating sublocalization is observed in the 2d case. 

Equation ( ^0| ) is supposed to be rigorous for sufficiently 
large r values, where the distribution of is exactly de- 
scribed by the log-normal distribution H\^{A, r). The decay 
of 5'Ar(r) can be approximated by Eq. (Ila) with an effec- 
tive localization exponent dif, which depends on r and N . In 
Fig. ^, we trace dq, as a function of r/A for several values of 
N using the parameters A(r) and a corresponding to the An- 
derson model in Id and 2d. We determined the values of 
by local fits of Eq. ( [T^ to Eq. (pO[). Surprisingly, the behavior 
is very different in Id and 2d. 

In Id, we find that the effective exponent d^ in the expo- 
nential ansatz (|l5|) is larger than one except for the typical 
average {N = 1) which we will discuss later. The localization 
is stronger than the usual exponential localization. Thus, this 
regime can be described by the term "superlocalization". Fur- 
thermore, the values of strongly depend on the number of 
configurations N . We observe a local maximum in the depen- 
dence of d^ on r for r w lOA to lOOA. While the height of the 
maximum increases with increasing N, indicating stronger 
superlocalization for an increasing number of configurations, 
the position of the maximum is shifted towards larger r val- 
ues. Only in the asymptotic limit for extremely large r values, 
r K, lO^A, the effective localization exponent d^, converges 
to one, yielding the theoretically predicted simple exponential 
behavior 

The typical average 5'i(r), describing the behavior of one 
typical eigenfunction, can be deduced by setting N = 1. As 
can be seen in Fig. ^ the exponent d^j is smaller than one 
for intermediate distances r from the localization center (for 
r < 60A). For larger r it is very close to, but slightly above 
one. Practically, this behavior can no longer be distinguished 



from simple exponential localization for r > 50A. 

In 2d, however, the localization behavior appears to be re- 
markably different. As can be seen in Fig. pi we observe 
hardly any A^-dependence of the effective localization expo- 
nent d*, although the absolute value of the average wave func- 
tion does slightly depend on N . The local maximum of is 
so weak, that it will be practically invisible. The most impor- 
tant difference to the localization behavior in Id, though, is 
the small value of d*. The effective localization exponent is 
smaller than one not only for intermediate r values, but even 
remains below one asymptotically up to at least r = 10^ A. 
Thus simple exponential behavior of the average wave func- 
tion can apparently not be reached for any practical system 
size. This is mainly due to the logarithmic dependence of the 
localization length A on r, that we have reported in Fig. ^ and 
inserted into Eq. ( pO[ ) to derive d^,. Consequently, this regime 
is not characterized by superlocalization as in Id. We rather 
find a stretched exponential decay, so that the term "sublo- 
calization" seems to be appropriate, because the localization 
behavior is weaker than simple exponential. 

Equation ( ^0| ) is not valid for small r values where the 
real distribution H{A^r) deviates from the log-normal form 
H\a{A,r). Thus, deviations occur for r approximately 
smaller than lOA. Further deviations occur for intermediate 
r and large N values when the tails of the distribution become 
more important. Then the lower integration limit A^i^{r, N) 
in Eq. (|l9|) tends to zero and the distribution of rare events of 
wave function amplitudes close to the maximum (corre- 
sponding to A= — 111 I V-" I close to zero) becomes important. 

Despite of this, the log-normal approximation H\a{A,r) 
can be used for obtaining a qualitative picture of 'l'jy(r) also 
at small r values, where Amiii(r, A) = in Eq. (|l8|). We 
introduce the cutoff at A = 0, since no amplitude ifju at dis- 
tance r > from the localization center no can be larger than 
iV'iio I ~ 1' ™d thus A = — In |-0„| > 0. It can be shown nu- 
merically, that the integration of Eq. ( (1^ ) yields an effective 
localization exponent much smaller than one, d^, « 0.5, for 
small r in the ansatz Eq. ([T5|). Thus, the decay is weaker than 
exponential in the small-r regime, and we expect a regime of 
sublocalization. With increasing A^, yl,„i„(r, A^) in Eq. ( p^ 
decreases reaching zero also for larger r values. Thus, the 
small-r regime with A,„in(r, A) — becomes dominant for 
averages over a large number of configurations. Since the 
lower integration limit becomes zero in Eq. ( p^ for the small- 
r regime, it is evident that ^ jv ('') does not depend on A^. 



V. NUMERICAL RESULTS FOR THE MEAN WAVE 
FUNCTION 

A. Averaging procedure 

Next we determine numerically the way the average ampli- 
tude values decay with increasing distance r from the 
localization centerx-JCba-ttecessary averaging procedure con- 
sists of three steps. c3l22l'E2l 
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FIG. 6.; The decrease of the typical average amplitudes for A'^ con- 
figurations, 'l'jv(r-), versus r for electronic wave functions in the 
Anderson model in Id for w = 1.0 and N = 1,16,256,4096, 
and 65536 (from the bottom to the top). For the plot, more than 
lO'' eigenfunctions have been calculated on chains of lengths L = 
30000. The straight lines with the slopes = 1 (dashed) and 
= 0.4 (dotted) are shown for comparison. In the inset, the ef- 
fective local exponents d* determined numerically from the slopes 
of the curves are shown versus r. The symbols correspond to the ef- 
fective numbers of configurations, N = 1 (circles), 16 (triangles), 
256 (squares), and 4096 (stars). 

1 . In the first step we average, for each eigenfunction ipn^ 
separately, the values of It/'rT'' | for all sites n at given distances 
r from the localization center. The resulting function 'ij}''^'^ (r) 
characterizes the spatial decrease of the v\h eigenfunction. 

2. In the second step, for obtaining the mean spatial de- 
crease of N eigenfunctions, we average over N con- 
figurations (eigenfunctions), 




FIG. 7.: The decrease of the typical average amplitudes for A*' 
configurations, 4'jv(r), versus r/\ for electronic wave functions 
in the Anderson model in Id for iV = 1, 16, 256, 4096, and 65536 
(from the bottom to the top) and five different strengths of disorder 
{w = 0.35, 0.5, 1.0, 2.0, and 4.0). For the total number of config- 
urations, the values of \{w), and the system sizes, see Fig. 2. The 
straight lines with the slopes d* = 1 (dashed), d* = 0.4 (dotted), 
and d* = 0.3 (dot-dashed) are shown for comparison. In the in- 
sets, the effective local exponents d* determined numerically from 
the slopes of the curves, are shown versus r/A. The symbols corre- 
spond to the effective numbers of configurations A'^ = 1 (circles), 
16 (triangles), 256 (squares), 4096 (stars), and 65536 (diamonds). 

The special case = 1 in the definition corresponds to the 
"typical average" over one configuration. Since = 
'ip('^)(r) in Eq. (pi]), the second step in the averaging proce- 
dure is dropped, and the values from each configuration are 
averaged logarithmically. 



Wr))^^-^^('')(r). (21) 

If the system is not selfaveraging (which we find is the case 
here), the resulting values for {ip{r)) jy will depend on the spe- 
cial set of N configurations considered and will fluctuate from 
set to set. 

3. To obtain the typical value of (V'(r))jv, we are led, in 
the third step, to the logarithmic average over many sets of N 
configurations, i. e_we. average ln(?/'(r)) jv over many sets of 
N eigenfunctions,E3 E1E3 

^N{r) = cxp(ln(V'(r))Ar). (22) 

In the following, we shall discuss exclusively this "typical av- 
erage over N configurations" defined by Eq. (p2|). 



B. Results in Id 



Figure^ shows 5'Ar(r) for the electronic wave functions in 
Id for several A^ values. The eigenfunctions are calculated fpt 
states in the band center (E « 0) by the iteration method,E3 
and the disorder parameter was chosen to be w = 1.0. As can 
be seen in the Figure, the actual values of 5'Ar(r) depend sig- 
nificantly on the number of configurations N included in the 
averaging procedure. The local slopes of the curves yield the 
effective exponent in our ansatz Eq. ( |l5[ ) which are shown 
in the inset of Fig. |[ One can see that the decay of "^Nir) is 
not simple exponential {d^ = 1), but rather two different lo- 
calization regimes with evidently different effective exponents 
can be distinguished. 
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FIG. 8.: The crossover distances rx{N)/X for electronic wave 
functions in the Anderson model in Id are shown versus A'^ in a 
semi-logarithmic plot. The values of rx (A'^) are determined as the 
lengths r, where the effective localization exponent of Vl/iv(r) 
intercepts with 1. The symbols correspond to the different disor- 
der strengths, w = 0.35 (diamonds), 0.5 (stars), 1.0 (squares), 
2.0 (triangles), and 4.0 (circles). The straight line rx{N)/X = 
3 logjQ A'' + 4 is shown for comparison. 

In the first localization regime, where ^'Ar(r) decays 
roughly by one order of magnitude, we find a stretched ex- 
ponential decay (sublocalization) with an effective exponent 
dij/ = 0.4 for all values of N, which is consistent with the re- 
sult obtained in the previous Section. As expected, this sublo- 
calization regime expands as the number of configurations A^ 
increases. 

In the second regime, for large r values, self-averaging fails 
and ^' jv (r) depends explicitly on A^. In the beginning of this 
regime, is considerably larger than 1, corresponding to a 
faster than exponential decay (superlocalization). Only for ex- 
tremely large r values the simple exponential decay (d^ = 1) 
is reached. This behavior agrees nicely with the analytical 
results obtained in the previous Section, as can be seen by 
comparing the inset of Fig. ^ with the analytical values of the 
effective localization exponent d^J, shown in Fig. || Note that 
very large r values cannot be reached in the numerical calcu- 
lations. Thus, the theoretical description, that becomes exact 
asymptotically for large r, is complementary to the numerics 
in this sense. 

The localization behavior is similar for other strengths of 
disorder w. This can be seen in Fig. 0, where the data for five 
different w has been included. Due to the scaling of the local- 
ization length, A — 105/u''^, the data can be scaled onto mas- 
ter curves for each number of configurations A^. For the sec- 
ond regime, where superlocalization occurs, we get an excel- 
lent data collapse. In the first regime, where sublocalization 
occurs, the data collapse is less accurate, because the local- 
ization exponent depends on the disorder strength w. We 
obtain dq, « 0.3 for small disorder (w = 0.5) and 0.4 
for larger disorder (w — 2.0). 

The distance characterizing the crossover from the first 
(r < Tx) to the second (r > ry) localization regime de- 
pends on the number of configurations A^. Since < 1 
in the sublocalization regime (for r < rx) and > 1 in 
the superlocalization regime (for r > rx), the values of rx 
can be determined as the lengths r where intercepts with 



unity. Figure ^ shows rx {N)/X versus log A^. It is obvious 
that rx (A^) increases approximately logarithmically with N, 

rx(^)/A«31ogio^ + 4. 



C. Results in 2d 



Figure ^ shows our results for the mean amplitudes jv (r) 
[defined in Eqs. ( pi] ) and (|2^] for the Anderson model in 
2d. We consider the same two disorder strengths and types 
of boundary conditions as in Figs. |^ and ^ and four values of 
A^. The Figure shows clearly that at all distances from the lo- 
calization center stretched exponentials (sublocalization) oc- 
cur, which is fully consistent with the theoretical results ob- 
tained in Section IV (see especially Fig. In addition, there 
is practically no dependence of the effective localization expo- 
nent, obtained from the local slopes in Fig. ^ on the number 
of configurations A^ taken into account in the averaging pro- 
cedure. The mean wave function amplitudes are just slightly 
shifted to larger values with increasing A^. 

In the small-r regime, which cannot exactly be described 
by the theoretical study in Section IV, we find strong sublo- 
calization characterized by « 0.33, similar to our results 
in Id where we also obtained values between 0.3 and 0.4 (see 
Figs. ^ and ||). This regime ranges over approximately two 
orders of magnitude in r from the site neighboring the center 
of localization (r — 1) to rx (A^), which is > 100 for large N 
in the two examples considered. In this first regime the mean 
wave function amplitudes decay to values of w 10^^. Thus, 
this regime is most relevant for most applications. 




FIG. 9.: The decrease of \l/jv (r) versus r for electronic wave func- 
tions in the Anderson model in 2d for TV = 1, 8, 64, and 512 (from 
the bottom to the top). For (a) w — 8.5 and HWBC were chosen, 
while w = 10.0 and PBC in (b). For each part of the figure, ~ 1000 
eigenfunctions with i5 ~ were calculated using the Lanczos al- 
gorithm on a 300 x 300 lattice. The straight dashed and dotted lines 
have the slopes dg< = 1 and = 0.33, respectively, and are shown 
for comparison. 
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The second regime for large distances r from the locaHza- 
tion center (r > r^) can be seen only in a quite limited r 
range, because of the limited system size and finite size ef- 
fects especially for PBC. The effective locahzation exponents 
dis, in this regime are definitely larger than in the first regime, 
but we still obtain values remaining strictly below one. Thus, 
in 2d, the second regime is also characterized by sublocahza- 
tion, in contrast to the superlocalization observed in Id. This 
numerical result is in very nice agreement with the prediction 
from the theoretical description in Section IV. 



VI. SUMMARY 

In order to discuss the localization and fluctuation behavior 
of localized eigenfunctions in the Anderson model we have 
studied the histogram amplitude distribution function H{A,r) 
which corresponds to the probability density to find the ampli- 
tude A = — In I V-" I at distance r from the localization center. 
We find that H{A, r) quite exactly follows a log-normal shape 
for Anderson wave functions in Id and 2d. For increasing dis- 
tance from the localization center, the log-normal fits become 
even better. Hence, the distributions can be characterized by 
two parameters, the localization length A and the fluctuation 
parameter a that is proportional to the width of the distribu- 
tions. While both parameters exactly follow the predictions 
of the standard single parameter scaling theory in Id, devia- 
tions occur in higher dimensions. The fluctuation parameter is 
smaller than the standard value of the single parameter scaling 
theory for the Anderson model in 2d, in contrast to a larger 
value for the percolation model. Furthermore, the localiza- 
tion length shows a logarithmic dependence on the distance 
from the localization center for the Anderson model in 2d, 
which gives an explanation for the contradicting values that 
have been published before and casts doubts on the existence 
of a finite asymptotic locahzation length. 

Using the log-normal ansatz for the amplitude distributions, 
the decay of the average eigenfunctions can be calculated ana- 
lytically. It is remarkable that by this simple approach, the es- 
sential features of the locahzation phenomenon, in Id: sublo- 
calization in the regime of small distances from the localiza- 
tion center, crossover to superlocalized behavior (that depends 
logarithmically on the number of configurations A^), and final 
approach to simple exponential behavior, and in 2d: sublocal- 
ization in all regimes, are reproduced. 
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